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Abstract 

10 — 1O 5 M black holes with dark matter spikes that formed in early minihalos and still 
exist in our Milky Way Galaxy today are examined in light of recent data from the Fermi 
Gamma-Ray Space Telescope (FGST). The dark matter spikes surrounding black holes in our 
Galaxy are sites of significant dark matter annihilation. We examine the signatures of annihi- 
lations into gamma-rays, e + /e~, and neutrinos. We find that some significant fraction of the 
point sources detected by FGST might be due to dark matter annihilation near black holes in 
our Galaxy. We obtain limits on the properties of dark matter annihilations in the spikes using 
the information in the FGST First Source Catalog as well as the diffuse gamma-ray flux mea- 
sured by FGST. We determine the maximum fraction of high redshift minihalos that could have 
hosted the formation of the first generation of stars and, subsequently, their black hole rem- 
nants. The strength of the limits depends on the choice of annihilation channel and black hole 
mass; limits are strongest for the heaviest black holes and annhilation to bb and W + W~ final 
states. The larger black holes considered in this paper may arise as the remnants of Dark Stars 
after the dark matter fuel is exhausted and thermonuclear burning runs its course; thus FGST 
observations may be used to constrain the properties of Dark Stars. Additionally, we comment 
on the excess positron flux found by PAMELA and its possible interpretation in terms of dark 
matter annihilation around these black hole spikes. 



1 Introduction 



The very first generation of stars, known as Population III. 1, likely formed from metal-free, molec- 
ular hydrogen-cooled gas at the center of dark matter minihalos of ~ 1O 6 M at z > 10 03-[6]|. 
Simulations indicate that they were quite massive, typically more than 100 M Q . If dark matter 
is made of Weakly Interacting Massive Particles (WIMPs), which are typically their own antipar- 
ticles, then the first phase of stellar evolution may be a Dark Star (DS) phase 0, during which 
the star is powered by dark matter (DM) annihilations prior to nuclear fusion. Indeed, if the DS 
phase persists for an extended period of time before fusion sets in, the star may grow to be cor- 
respondingly more massive: as long as dark matter annihilation powers the star, it remains cool 
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enough to continue to accrete matter and grow ever larger. Initial work showed that they can grow 
to be ~ 1OOOM jH, while more recent studies, taking into account triaxial halos with a variety of 
particle orbits, demonstrate that they may even become super-massive > 1O 5 M 10. Subsequent 
to the initial paper on Dark Stars [7], other authors have explored the repercussions of DM heating 
in the first stars, including [fTQ] - l24| . 

When the DM fuel inside a DS runs out, it collapses and heats up to become a standard fusion- 
powered star. If some Population III. 1 stars experience a DS phase, the Initial Mass Function 
(IMF) for fusion-powered stars may be quite different than what might otherwise be expected Il25ll . 
Specifically, the IMF would be determined by the length of the DS phase, i.e. by how long it 
took each DS to exhaust its DM fuel supply. The "dark power" may have lasted different amounts 
of time in different DSs: millions to billions of years for some stars or, in extreme cases, even 
until today. Stars in the mass range 140 — 26OM would have ended their lives as pair instability 
supernovae, leaving no remnants H26U27H . but stars outside this mass range would collapse to black 
holes. It is the latter case that we consider in this paper. 

Our work, though motivated by the DS scenario, applies generally to black holes at the centers 
of minihalos, regardless of their origin. Many black holes of mass 10 — 1O 5 M that formed at 
the centers of minihalos survive in the universe today. Assuming some fraction of high redshift 
minihalos hosted Population III. 1 star formation, one can estimate the distribution of their remnant 
black holes today. Each remnant black hole will be surrounded by a region of enhanced dark 
matter density, which we refer to as a "spike." Although some of the original minihalos would 
have merged with other DM halos, resulting in disruption, one can still follow the evolution of the 
black holes and their DM spikes in simulations. Here, we use the Via Lactea-II simulation to track 
black hole spikes from the redshift of their formation to z = 0. In this way, we estimate the black 
hole population in a galaxy like our Milky Way today. 

When the baryonic gas in a minihalo collapses to form a structure, be it a star or black hole, 
dark matter is dragged in, creating a DM spike around the central object. The star or black hole 
provides a gravitational potential that causes the surrounding DM to be pulled inward, into and 
around the object. The resulting enhanced DM density may be computed via adiabatic contraction, 
assuming that the DM orbital timescale is shorter than the timescale of the changing gravitational 
potential. As the star reaches the end of its lifetime, collapsing to become a black hole, the DM 
spike remains in place. The enhanced DM density in these spikes leads to an enhanced rate of DM 
annihilation, which scales as the square of the DM density. Unstable DM annihilation products 
undergo hadronization and decay, with typical final products being gamma-rays, e + /e~ pairs, and 
neutrinos, each of which may provide a significant signal at a variety of ground- and/or space-based 
detectors. 

In this paper, we explore the detection prospects for gamma rays produced in dark matter 
annihilations in the DM spikes surrounding black holes for a range of star formation scenarios, 
black hole masses, and dark matter annihilation modes. The current data from the Fermi Gamma- 
Ray Space Telescope (FGST) are used in a two-pronged approach to constrain the number of black 
holes in the Milky Way halo and, consequently, the number of DSs that could have formed at early 
times. We use the FGST First Source Catalog ||28l to find the minimal distance to the nearest DM 
spike such that it is not brighter than the brightest source observed by FGST. From the predicted 
distribution of such spikes in the Milky Way halo, we extract a limit on the fraction of minihalos 
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in the early universe to host a black hole (and survive as a DM spike in our Galactic halo today). 
We also use the FGST measurement of the diffuse gamma-ray background [29 ] to constrain the 
population of DM spikes contributing to the diffuse flux today, thereby setting a second limit on 
the fraction of minihalos in the early universe to become DM spikes in our Galactic halo. 

Several groups have previously studied signatures from DM annihilation in DM overdensities, 
or spikes, around black holes. Gondolo and Silk ll3~0ll first pointed out the possibility of a DM spike 
around the ~ 1O 6 M Super-Massive Black Hole (SMBH) at the center of the Galaxy; they coined 
the terminology "spike" to differentiate the gravitationally contracted DM due to the existence of 
the black hole from the more standard possibility of cusps at the centers of DM halos in cold dark 
matter scenarios of structure formation. The enhancement in the dark matter density in the spike 
due to the SMBH at the center of the Milky Way was later shown to be likely reduced by a variety 
of effects including mergers, formation of the initial black hole away from the Galactic Center, 
and gravitational scattering with stars [[3TI434I . Zhao and Silk 11351 suggested DM spikes around 
Intermediate Mass Black Holes, remnants of Pop. Ill stars, similar to the work in this paper. Shortly 
thereafter, Bertone, Zentner, and Silk (BZS) (35|, with follow-up work in 1371. Q8l and 11391, 
studied formation and evolution histories for black holes and proposed looking for signatures of 
DM annihilation in the spikes around them in the Galaxy today. They chose two scenarios to 
examine: (i) 100 M black holes left over as remants from the first stars, and (ii) 1O 5 M black 
holes that formed in larger 1O 7 M halos due to accretion disks. They too estimated the distribution 
of unmerged spikes that reside in the Milky Way today for each scenario. Our work, although 
similar to theirs, differs quantitatively in the following respects. First, we identify star-forming DM 
halos in a different way, as described in Section[2] Second, since we are motivated by a potentially 
long-lived DS phase for Pop. III. 1 stars, we consider a variety of black hole masses between 10 
and 1O 5 M . Third, we use recently released data from FGST to constrain the current population of 
black holes and, consequently, the initial population of their progenitors. In particular, we consider 
data from both the FGST First Source Catalog [|28l and the measurement of the diffuse gamma-ray 
flux 11291 . 

In Section [2] we discuss the remnant black hole population in the Milky Way halo, and in 
Section [3J we discuss the DM spikes surrounding those black holes. In Section |4] we present the 
signal from dark matter annihilation in DM spikes in the Milky Way halo, that due to spikes that 
appear as point sources and the expected diffuse gamma-ray flux. We then use FGST point source 
data and the measurement of the diffuse gamma-ray background to place constraints on the star 
formation history in a variety of dark matter annihilation models, as presented in Section \5\ In 
Section [6] we briefly address signatures in positrons and neutrinos and discuss prospects for future 
work. Finally, in Section [7] we present our conclusions. 

2 Remnant Black Hole Population in the Milky Way Halo 

In order to estimate the observable consequences of annihilation in the DM spikes around black 
holes, we need to estimate the number of black holes in the Milky Way halo today that are the 
remnants of the first stars. Ideally, we would like to know the Initial Mass Function and Star 
Formation Rate of Pop. III. 1 stars; however, what guidance we can get from simulations is limited 
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by our understanding of feedback processes in the early universe. During the DS phase, since the 
stars are extremely cool there is very little ionizing radiation and therefore very little feedback. 
However, as soon as the first fusion-powered stars form, they produce ionizing photons which 
influence and hinder the formation of stars in neighboring regions. The effects of feedback on star 
and/or protogalaxy formation have been studied in Refs. H0] - l45l . among others. Because of the 
many uncertain aspects of theoretical studies of Population III. 1 star formation, in this paper we 
follow the method discussed in the remainder of this section. 

We use the Via Lactea-II simulation to estimate the number and mass distribution of DM mini- 
halos as a function of redshift. The minihalos that hosted Pop. III. 1 stars must have had masses 
in a very particular range at the time of star formation. Pop. III. 1 stars could only form when a 
cooling mechanism for the collapse of the baryonic cloud arose. The first accessible mechanism to 
cool the gas was via excitations of molecular hydrogen, the fraction of which present in a minihalo 
is related to the temperature, which in turn can be written in terms of the mass and redshift of 
the minihalo. We use the parametrization of Ref. [|46ll for the minimum halo mass in which star 
formation could occur: 



Because of the hierarchical nature of structure formation, there are far more smaller halos than 
larger ones. Our results are therefore not sensitive to the maximum halo mass for Pop. III. 1 star 
formation, which we take to be = 10 7 M Q . See also Iflli 

Next we assume that some fraction, f^ s , of viable minihalos hosted a Pop III. 1 star (e.g. a 
DS). We also assume that each star ended its life in collapse to a black hole. Neglecting, for now, 
black hole mergers, the comoving number density of black holes as a function of redshift is then 



where Nh a io{z) is the comoving number density of minihalos in which Population III. 1 star forma- 
tion was possible. 

The duration of the DS phase is highly variable, and DSs may live vastly different amounts of 
time before running out of DM fuel, becoming fusion-powered, and finally collapsing to a black 
hole. However, the lifetime of the star is not important as long as it has reached the black hole 
stage today. We also note that we do not consider further accretion of mass onto the black hole; 
this is the most conservative assumption regarding potential detectability and constraints on DS 
scenarios. 

At some time between the beginning of Population III. 1 star formation and the end of reion- 
ization at z ~ 6, hydrogen deuteride (HD) cooling would have become possible in minihalos, at 
which point massive Population III. 1 (and DS) formation would have given way to less-massive 
Population III. 2 star formation I1421 . however there are few constraints on when this transition 
occured. Motivated by recent discussions of the effect of feedback from the very first stars on 
subsequent star formation and reionization II4T14451I471 . we consider three scenarios for the ter- 
mination of Population III. 1 star formation at redshift Zf. These scenarios are hereafter noted as 
Early, Intermediate, and Late, following Ref. [42J, with Zf ps 23, 15, and 11, respectively, as shown 
in Tabled] For each case, we assume that Pop. III. 1 star formation was only possible in minihalos 
with masses in the range discussed above and prior to the termination redshift. 




(1) 



N BH (z) = f° DS N halo (z), 



(2) 
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Label 


z f 


Representation 


Early 


23 


green 


Intermediate 


15 


red 


Late 


11 


blue 



Table 1: Three scenarios for the redshift at which Population III. 1 star formation ceases, Zf. In the 
third column we list the color used to represent each case in the following figures. 



We determine the z = distribution of DM spikes throughout the Galactic halo from the Via 
Lactea II (VL-II) cosmological N-body simulation [48]. With a particle mass of 4.1 x 1O 3 M , 
VL-II is the first simulation of the Galactic DM halo that is able to directly resolve the small 
progenitors (minihalos) which are the formation sites of Population III. 1 stars. 

As discussed above, since the truncation redshift for Pop. III. 1 star formation is poorly con- 
strained, we consider three different scenarios: Early (zj = 23.1), Intermediate (zj = 14.8), and 
Late (zj = 11.2). At these redshifts we identify all minihalos between the minimum mass in Eq.Q] 
and a maximum mass of 1O 7 M Q and assign one Population III. 1 star to its most bound particle 
(tracer particle). We assume that this Population III. 1 star will then form a black hole surrounded 
by a DM spike, as described in Section [3] and that both the black hole and the spike will survive 
until z = 0. The z = position of the tracer particles directly relates to the distribution of DM 
spikes in this scenario. The simulation naturally includes the effects of dynamical friction and 
tidal stripping as the minihalos fall into larger host halos. Eventually many of these minihalos fall 
below 10 5 M Q and are no longer resolved in VL-II. From this point on the remaining minihalo with 
its black hole and DM spike is represented as a point mass in the simulation, i.e. the dynamical 
friction of the stripped minihalo against the Galactic DM halo is neglected, but this effect is very 
small for objects below 1O 5 M even near the Galactic center [50]. 

We mention here two caveats: First, VL-II does not include the formation of the Galaxy which 
could contract the DM halos and therefore also the spike distribution. However, contraction may 
be much smaller than the classical adiabatic models assume, if it exists at all (e.g. Ref. ||49l argues 
for DM halo expansion during galaxy formation). Second, when two minihalos (each containing 
a black hole) merge, it is possible that the black holes form a close binary which would destroy 
the DM spikes. Our model does not account for that. However, we estimate, as described in detail 
below, that mergers may change the number of DM spikes in our Galaxy today by at most a factor 
of 2 for the highest black hole masses considered, and very little for low mass black holes. 

The limited time and spatial resolution of VL-II and the lack of massive black hole particles in 
this DM-only simulation make it impossible to identify black hole mergers directly. To estimate 
the fraction of DM spikes which could have been destroyed in mergers we simply identify for each 
spike its nearest companion black hole in 24 snapshots between 2 = 11 and 2 = 00. For simplicity 
we assume here that the spikes orbit around each other on a circular orbit within an isothermal 

1 Using a larger maximum mass of 10 s M© gives virtually the same results, a consequence of the steepness of the 
halo mass functions for these redshifts and mass ranges. 

2 A list of the exact reshifts is available at www.physik.uzh. ch/^diemand/vl. 
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sphere, in this case the dynamical friction time is simply 11501 

t 19Gyr/ Ar \ 2 Av 1O 5 M 

In A V°- 5k P c / 20 kms" 1 M BH +DMspik e ' 

where Ar and Av are the distance and relative velocities of the two neighbouring black holes. 
The Coulomb logarithm is set to In A = 6. Equation [3] is evaluated for all DM spikes at each 
of the 24 time steps between z = 11 and z = and all DM spikes which might have merged 
by today according to their £f ric are added up. The fraction of DM spikes potentially lost through 
mergers is significant only for the most massive black holes and for f DS close to one: For the 
Late scenario, with f^ s = 1 and M B H+DMs P ikc = 1O 5 M the fraction of possibly lost spikes is 
f merged = 0.49. Lowering the mass to 1O 4 M results in f merged = 0.24, and going down to 1O 3 M 
gives f merged = 0.076. These estimates illustrate that these mergers would result in disruption of, 
at most, fewer than half of all Milky Way black holes for the largest black hole masses considered 
here, and much smaller fractions for lower black hole masses. Since it is not possible to directly 
identify black hole mergers in VL-II, for the following analysis we define fns to be the fraction of 
surviving black holes with surrounding DM spikes, such that it is related to the initial fraction of 
star-forming minihalos, fp S , as 

fDS — /ds(1 — f merged), (4) 

with the fraction of DM spikes destroyed in mergers, f merged, itself a function of fp S . For small 
f% s , mergers become very rare and fos — > Ids- 

In Fig. Q] we show the number densities of DM spikes inside the Milky Way halo as functions 
of Galactic radius for Early, Intermediate, and Late Zf. In the Early case, Pop. III. 1 star formation 
terminates at z m 23, so there were the fewest stars, and therefore the fewest black holes and 
surviving density spikes today. In the Intermediate and Late cases, Population III. 1 star formation 
turns off at redshifts of roughly 15 and 11, respectively. For comparison, the total dark matter 
density profile at z = in VL-II is also shown; although the normalization of these points is 
arbitrary, it is useful to illustrate that the total DM profile is more extended than the distribution of 
black holes with DM spikes. 

We can contrast our approach for finding the relevant DM minihalos to that of Ref. If3~6l . At 
z = 18, they populated halos that constituted 3a peaks in the smoothed primordial density field 
with seed black holes of initial mass 100M Q . Using an analytical model of halo evolution, they 
simulated 200 statistical realizations of the growth of a Milky Way- sized halo. A fit to their re- 
sulting distribution of DM spikes is shown as the grey curve in Fig. [TJ Instead, we use a single 
iteration of the VL-II cosmological simulation and follow potential Population III. 1 star- forming 
minihalos within the halo mass range and redshift ranges discussed above. 

Using their approach, Ref. Il36ll found that 1027 ± 84 unmerged EVIBHs and surrounding DM 
spikes are expected to exist in the Milky Way halo today. Our simulations yield 409, 7983, and 
12416 DM spikes in the Early, Intermediate, and Late scenarios, respectively, assuming fr>s = 1- 
However, not every viable minihalo must have hosted a Pop. III. 1 star, in which case the total 
number of DM spikes in our Galactic halo would be roughly proportional to fp S . Again, for small 
fp S , the fraction of DM spikes destroyed in mergers becomes negligible, and f DS w f^ s . 
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Figure 1 : The number density of black hole spikes in the Milky Way as a function of Galactic radius 
for star formation models with Early (green), Intermediate (red), and Late (blue) Zf as described 
in the text and for fos — 1- Curves have been obtained using the VL-II N-body simulation as 
described in the text. The black points illustrate the total dark matter density profile at z = in 
VL-II. Also shown as a solid grey curve is the analytical fit found in Ref. I[36\\38\l . Our simulations 
show 409, 7983, and 12416 DM spikes in the Milky Way for the Early, Intermediate, and Late 
scenarios, respectively, assuming fos — 1- 



3 Dark Matter Density Spikes 

The density profile of a dark matter spike surrounding a black hole is determined by adiabatic 
contraction of the dark matter halo around the central mass. For concreteness, as our starting point 
we take the Navarro, Frenk, and White (NFW) profiles for both the baryons (15% of the mass) and 
dark matter (85% of the mass): 

r/r s (l + r/r s ) z 



where p is me scale density and r s is the scale radius H5TH . The scale density, p , can be re- 
expressed in terms of the critical density of the universe at a given redshift, p c (z), via 

,200 C 3 

pQ - pc[z) ~in{i + o-c/{c + iy (6) 

where we take the value of the concentration parameter to be C = r vir /r s = 3.5, and r vir is the 
virial radius of the halo. We find that decreasing the concentration parameter to C = 2 results in a 
~ 30% decrease in the luminosity of each spike. Sensitivity of DS properties to the concentration 
parameter is discussed in Ref. Il52l . 
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We assume a flat ACDM universe with current matter density Vt m = 0.24 and dark energy 
density £7a = 0.76. 

We allow a point mass to grow at the center of the halo, mimicking the formation of a star 
and the subsequent remnant black hole. To model the response of the DM to this density growth 
we have used the simple Blumenthal et al. prescription for adiabatic contraction [53]. It has been 
shown that this method obtains a DM density profile that is accurate to within a factor of 2 [fT2l . 
The result of the adiabatic contraction is a roughly power law density profile. 

In Fig.|2]we show the contracted halo profiles today for DM spikes due to black holes of various 
masses for the case where the central object formed at z = 15. We note that the power law portion 
of the profile is independent of WIMP mass. In the central regions, closest to the black hole, some 
of the DM has annihilated away in the time since the formation of the central mass. We follow 
BZS who found an upper limit to the DM density 



where tsH is the lifetime of the central mass, roughly 1.3 x 10 years for a star that formed at 
z = 15. Thus we take the density profile in the inner region to be a flat plateau with this density out 
to the radius r*, defined by p(r*) = p max , beyond which the density profile follows that expected 
from adiabatic contraction of the minihalo. We note, however, that a proper treatment of this inner 
region is more complicated: Because dark matter particles in triaxial halos may follow box or 
chaotic orbits ll54l . some of the DM particles that are expected to have annihilated with each other 
in the center of the spike would not have been on orbits that would have returned them to the center 
anyway. Thus, for triaxial halos, it is harder to deplete the density in the central region. The use of 
Eq. [7] may therefore underestimate the amount of DM remaining in the center, but on these long 
timescales it is an acceptable estimate. 

Fig. |3] illustrates the dependence of the DM density profile on the redshift of formation of the 
central DS and subsequent black hole. The right panel reveals that there is nearly a factor of two 
difference between the density profiles for black holes that form at z = 10 and z = 20 at the radii 
shown, which correspond to regions that contribute significantly to the total DM annihilation rate 
(we note that only a fraction of a percent of the total annihilations in each spike occur at r > 1 pc). 
The reason for this difference is the fact that at higher redshifts, densities are higher as (1 + z) 3 
and the halos have smaller radii (n.b. we keep C fixed throughout). Thus, we expect differences 
in the flux of annihilation products from DM spikes created at different redshifts, all other factors 
being equal. 

4 Gamma Ray Signal from Dark Matter Annihilations 

For a Majorana dark matter particle with mass m x and annihilation cross section times velocity 
(crv), the rate of WIMP annihilations in a DM spike is 



m 



Pmax 



(<Tv)t B H 



(7) 
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Figure 2: Density profiles for contracted dark matter halos surrounding black holes of mass 1OM 
(magenta), 1O 2 M (blue), 1O 3 M (red), 1O 4 M (green), and 1O 5 M (orange), from bottom to top, 
assuming m x = 100 GeV and that the central black hole in each case formed at z = 15 in a halo 
of mass 1O 6 M . Note that 1 pc = 3.1 x 10 18 cm. 

with r min and r max defining the volume of the DM spike in which annihilations occur. At a 
minimum, r min should be equal to the Schwarzschild radius of the black hole, r Sch = 2Gitibh 
for a black hole of mass m BH and Newton's constant, G. We take r min = Ar Sc hi though our 
results are not sensitive to this choice for any r min < 10 14 cm, corresponding to ~ 3 x lO^ 8 ^ 
Schwarzschild radii for central black hole masses of 1O*M . 

We choose as a benchmark scenario (av) = 3 x 10~ 26 cm 3 s _1 , in agreement with the measured 
dark matter abundance today for thermal WIMP dark matter, and consider several WIMP candi- 
dates, defined by mass and annihilation channel. Calculations are performed for WIMP masses 
of 100, 200, 500, 1000, and 2000 GeV and Standard Model final states 66, W + W~, r + r~, and 
The resulting spectrum of photons from annihilation to final state /, dNf/dE, is computed 
with PYTHIA [|55l . For \X ~^ A t+ A i ~» the photon spectrum comes from final state radiation and is 
given by 

dNn +li - _ a fx 2 -2x + 2 

dx 7T \ X 

where x = E^jm x , the center-of-mass energy squared is s = 4m£ , and a ~ 1/137 [|56l . We note 
that WIMP candidates such as neutralinos typically annihilate to a variety of final states, /, with 
branching fractions Bf such that the total annihilation cross section (av) can be expressed as the 
sum over final states 

m=x;w>/= w>E s /« (10) 
/ / 

with the final sum evaluating to unity. Consequently, the rate of annihilations to final state / in a 
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Figure 3: Density profiles for contracted dark matter halos surrounding black holes of mass 1O 2 M 
for black hole formation at z = 10 (blue, dashed), z = 15 (red, solid), and z = 20 (green, dotted). 
The initial halo mass is taken to be 10 6 M & for each of the three cases, and C = 3.5. The right panel 
zooms in to a range of radii that contribute significantly to the total DM annihilation rate in the 
spike to show more clearly the differences between the three cases. Note that 1 pc = 3.1 x 10 18 cm. 

DM spike may be expressed as 



The intrinsic photon luminosity from dark matter annihilations in any dark matter spike is then 



with T j given by Equations [8] and [TTJ This quantity is plotted as a function of WIMP mass in 
Fig. |4] assuming Bf — 1 for each of the final states bb (blue), W + W~ (orange), t + t~ (magenta), 
and (green) and for DM spikes surrounding black holes of mass 1OM (thin dotted curves), 

IOOMq (thick solid curves), and 1000M o (thin dashed curves). We see that increasing the mass 
of the central black hole by an order of magnitude increases the luminosity by nearly an order of 
magnitude for all WIMP masses and final states. Also, heavier WIMPs result in significantly lower 
luminosity, as there are fewer of them in each DM spike and therefore the annihilation rate is lower. 
We note also the similarity of the luminosities for final states bb and W + W~, a consequence of the 
similarity of the photon spectra dN/dE in these two cases. 

For simplicity, in the following sections, we focus on 8 combinations of WIMP mass and final 
state, which we label with one particle of the final state and the WIMP mass, as shown in Table [2] 
Furthermore, we assume that dark matter has only one dominant anihilation mode, and consider 
only Bf = 1 for each final state. If dark matter annihilations result in more than one final state, 
as for standard neutralino or Kaluza-Klein dark matter, our results must be scaled appropriately. 
Given a WIMP annihilation model, the distribution of DM spikes in the Milky Way halo (as shown 
in Fig.Q]), and the characteristic density profile of each DM spike (for example, as shown in Figs. [2] 
and©, we calculate the expected gamma-ray signal. In Section l4Tl we discuss FGST point sources, 
and in Section H31 we discuss the diffuse gamma-ray flux. 



r/ = BjT. 



(11) 




(12) 
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Figure 4: The intrinsic luminosity of photons with energy above 1 GeV from a single DM spike as 
a function of WIMP mass in the Intermediate Zf scenario. The thick curves are for a central black 
hole of 10 2 Mq and the thinner dotted and dashed curves are for central black holes oflOM & and 
10 3 M©. From top to bottom, the blue, orange, magenta, and green contours represent final states 
bb, W + W~, t + t~, and /i + pT, respectively. 
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1000 


bb 


rlT 


1000 
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W100 


100 
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//100 


100 




WIT 


1000 


w+w~ 


/ilT 


1000 





Table 2: example WIMP annihilation models 



4.1 Point Sources 

The differential flux of neutral particles from annihilations to final state / in a DM spike with 
radius r max located some distance D from our Solar System is given by 

d$ f T f dN f 

— - = — (13) 

dE AnD 2 dE' 

for D 3> r max . If D is not large compared to r max , however, an integral must be performed along 
the line-of- sight, s, and over the solid angle of interest on the sky, defined by the polar angle, 9: 



d$ f _ B f (av) dN 



dE 2m 2 x dE 



j d0 27rsin0 / dsp 2 DM (r), (14) 
Jo Jo 
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where r = \J D 2 + s 2 — 2sD cos s ma:r > D + r ma;r , and 9 max defines the solid angle to which 
the detector is sensitive as Af2 = 2n(l — cos 9 max ). 

FGST has compiled a catalog of point sources identified at > 4cr significance within the first 
year of data collection 11281 . If a single DM spike is a sufficiently bright and compact source of 
gamma-rays, it may have been identified as a point source and recorded in the FGST First Source 
Catalog, which contains 1451 point sources, including 630 that are not associated with sources in 
other astronomical catalogs. Here we consider only point sources that are more than 10° away 
from the Galactic plane (that is, \b\ > 10°), as it is expected that the region close to the plane will 
contain the majority of baryonic gamma-ray sources (pulsars, supernova remnants, X-ray binaries, 
etc.). Of the unassociated point sources in the FGST First Source Catalog, 368 have been detected 
with greater than 5cr significance and are more than 10° away from the Galactic plane. 

Ref. [57] identifies a set of criteria for objects which should be identified as point sources by 
FGST in the first year of observation, which we adopt here, as wel§ Specifically, > 50 events 
per year must be observed by FGST and > 95% of the events must come from within a cone of 
half-angle 2° centered on the source. A bound on our black hole sources can be placed as follows: 
spikes bright enough to be identified as point sources must not be brighter than the brightest source 
in the the FGST First Source Catalog, which has a flux of photons with 100 MeV < £ 7 < 100 
GeV of 1.25 x 10~ 6 cm~ 2 s _1 . Requiring that the flux not exceed this value establishes a minimal 
distance, D^f n , beyond which the spike must be located in order not to be brighter than any source 
in the FGST catalog. 

The DM spikes from adiabatically contracted minihalos around black holes have quite steep 
density profiles, so in all models considered we find that all sources located at distances > D^f n 
are highly localized such that the 2° requirement is more than satisfied. We would like to empha- 
size that the dark matter spikes in our paper are indeed point-like objects, with the dark matter 
annihilation signal coming from the spike rather than from the extended minihalo around it. Of all 
the scenarios considered, the closest a spike can be to our solar system is 0.01 kpc (see Fig. 5). 
Even in this case, > 99% of the signal comes from within a cone of half-angle 0.01° (that is, most 
of the annihilations occur within ~ .002 pc of the central black hole). Since the 95% confidence 
level Fermi containment angle is at least an order of magnitude larger than this at all energies ||58l , 
this nearest spike is decidedly a point-like object. 

Before proceeding, we would like to point out that the brightest FGST point source is actually 
identified as associated with the Vela pulsar. The brightest source that is unassociated with known 
sources in other wavelengths has a flux of 5.78 x 10~ 8 cm _2 s _1 . The factor of ~ 22 difference 
between the flux from the Vela source and that from the brightest unassociated source means that 
had we chosen to base D^ n on the unassociated source we would find each D^ n to be a factor of 
\/22 = 4.7 larger. We will return to this point in the following analysis. 

In Fig. |5l we display in the left panel the minimal distance, D^ n , at which a point source must 
be located in order not to exceed the largest flux from any point source measured by FGST for the 
eight example cases in Tabled Also plotted in Fig. [5] is a horizontal line indicating the distance 
from our Solar System to the Galactic center. Note that for 100 GeV WIMPs annihilating to bb or 

3 In order to identify point sources in the first year data, the Fermi Collaboration have used a complicated diffuse 
emmission model and likelihood Test Statistic as detailed in Ref. |28l . The sensitivity of our results to our choice of 
simplistic criteria is discussed in Secs. l4.2l andl5l 
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Figure 5: In the left panel, we display the minimal distance from our Solar System of a single spike 
such that it does not exceed the flux of the brightest source in the FGST First Source Catalog as 
a function of central black hole mass for Intermediate Zf. From top to bottom, the contours are 
for models blOO (blue dashed), W100 (orange dashed), rlOO (magenta dashed), blT (blue solid), 
WIT (orange solid), yulOO (green dashed), rlT (magenta), and /xlT (green solid). In the right 
panel, we show the maximal distance from our Solar System of a single spike such that it would 
appear as a > 5cr point source to Fermi for the same cases. We assume that Bj = 1 for each of 
the annihilation modes considered. If the total annihilation cross section is the sum of branching 
fractions to different Standard Model final states, D^f n and D^ ax would be scaled by ^/B~f. In 
each panel, the horizontal black line indicates the distance to the Galactic center. 



W + W~, if ttibh ?t 2 x 10 4 M Q the nearest spike to our Solar System must be at least as distant as 
the Galactic center. From the spike distributions in Fig. Q3 we see that this is extremely unlikely. 
In fact, within only 5 kpc of our Solar System, we expect there to be ~ 7 DM spikes for Early 
termination of star formation, and 94 DM spikes in both the Intermediate and Late cases, assuming 
fas = 1- Since the number of spikes in any volume scales with fos, we find that fos can be 
quite constrained for very large central black holes. However, for most of the DM models shown, 
if the central black hole is < 1O 3 M , spikes may be located within 1 kpc of our Solar System. No 
matter which of the three choices for z/ we adopt, we expect that less than one DM spike would 
be located within 1 kpc of our Solar System, even for f DS = 1. Again, had we chosen to base our 
analysis on the brightest unassociated point source, all curves in the left panel of Fig. [5] would be 
shifted up by a factor of ~ 4.7. 

In the right panel of Fig. [5l we show the maximal distance, D^ x , at which a single DM 
spike would be identified by FGST as a > 5a point source. These distances were determined by 
requiring 50 events per year in FGST from each spike. For 100 GeV WEVIPs annihilating to bb 
or W + W~ around large black holes, as in the upper right portion of the right panel of Fig. [51 all 
DM spikes in the Milky Way halo would be identified as point sources by FGST. In the opposite 
extreme, for annihilation of 1 TeV WIMPs to t + t~ or in the DM spike surrounding a 10 

Mq black hole, any spike further than ~ 1 kpc would not be bright enough to be identified as a 
point source. As mentioned above, it is not guaranteed that there are any DM spikes within 1 kpc 
of our Solar System; thus no DM spikes may be bright enough to appear as FGST point sources. 
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The results shown in Fig.[5jare for Intermediate Zf (star formation ending at Zf ~ 15). In the Early 
and Late Zf scenarios there is a ~ 15% decrease and a ~ 13% increase in D^ s ax , respectively, and 
similar changes in D^f n . 

In Table[3]we list the number of DM spikes that would have been identified as 5a point sources 
by FGST in the first year of operation in each of the WIMP annihilation scenarios in Table[2]for five 
black hole masses and Early, Intermediate and Late Zf. The number of such sources with Galactic 
latitude \b\ > 10° is shown in parentheses next to each entry. There are 368 unassociated 5a point 
sources at |6| > 10° in the FGST First Source Catalog Il28l 15711 . One can see that many models 
studied here greatly overproduce FGST point sources if fns = L but recall that the number of 
point sources scales with fns- We note also that not all sources with the requisite brightness have 
been detected by FGST, for example due to local background fluctuations or spectral shapes that 
did not pass all analysis cuts. Even with these considerations, it is clear from Table [3] that many 
conclusions regarding f DS are highly model-dependent. One can say with certainty, however, 
that f DS is most constrained for large central black holes in the Intermediate and Late Zf cases, 
and possibly entirely unconstrained for smaller central black holes and if Pop. III. 1 star formation 
terminated at Zf > 20. 

Additionally, we remind the reader that D^f n has been determined based on the brightness of 
the brightest source observed by FGST, rather than by the brightness of the brightest unassociated 
point source, thereby over-estimating the number of nearby point sources in Tabled If we assume 
that all gamma-rays from associated sources are in fact due to the objects with which they are 
associated (and not DM spikes), then any DM spike must not be brighter than the brightest unas- 
sociated source. D^ n would therefore increase by a factor of 4.7 and consequently the number of 
point sources between D^f n and D^ x would decrease. 

Furthermore, as we have defined D^ x based on the assumption that objects that generate 50 
events per year in FGST would be identifiable as point sources, if the number of required events to 
pick out a point source is actually larger than this, then the numbers of point sources in Table[3]are, 
again, too large. In fact, the diffuse gamma-ray background is not entirely uniform, and the Fermi 
Collaboration's ability to identify point sources in the data depends on complex modelling of the 
diffuse gamma-ray sky and the location of the sources in the sky ESI . Since the number of spikes 
that would appear as point sources depends on the flux required for each spike to be identified as 
a point source, and since this flux depends on the region of the sky in which the spike is located, 
we do not use our estimates of the numbers of spikes in our Galactic halo to constrain f DS . We 
encourage the reader to interpret the numbers in Table [3] merely as reasonable estimates, with the 
actual number of point sources likely being somewhat smaller. 

4.2 Diffuse Flux 

Many of the DM spikes in the Milky Way halo may be too faint to be identified as point sources 
by FGST. Indeed, even DM spikes that have large enough fluxes to be identified as point sources 
may have been missed as a result of selection effects; for example, if the local background near 
the spike is large ||59il . Those DM spikes that cannot be identified as point sources may contribute 
to the diffuse flux. Based on the point source criteria in Sec. 14. II and the simulated distribution of 
Pop. III.l/DS remnant spikes in the Milky Way halo from VL-II, we calculate the contribution to 
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195 (70) 


1117(649) 


557 (387) 


~ 2 (~ 2) 


21 (17) 


14(11) 


10 2 


304(151) 


3247 (2263) 


2935 (2186) 


147 (45) 


586 (372) 


281 (215) 


10 3 


380 (213) 


5715 (4283) 


6754 (5305) 


284(135) 


2788 (1895) 


2340 (1708) 


10 4 


381 (217) 


7237 (5548) 


10608 (8486) 


372 (207) 


5213 (3870) 


5866 (4575) 


10 5 


158 (128) 


5918 (4831) 


10946 (8946) 


392 (224) 


7069 (5402) 


9998 (7980) 
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Late 


Early 
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176 (58) 


868 (504) 


392 (289) 


~ 2 (~ 2) 


24 (20) 


16(13) 


10 2 


294 (144) 


3011 (2073) 


2618 (1930) 


159 (50) 


692 (420) 


321 (242) 


10 3 


377 (211) 


5461 (4072) 


6279 (4915) 


288 (139) 


2881 (1969) 


2457 (1801) 


10 4 


388 (222) 


7152 (5473) 


10318 (8242) 


374 (209) 


5320 (3957) 


6043 (4718) 
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169(136) 


6140 (4998) 
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391 (223) 


7107 (5434) 


10144 (8099) 
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36 (30) 
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< 1 (< 1) 


< 1 (< 1) 
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211 (81) 


1360 (808) 


768 (518) 


~ 1 (~ 1) 


13 (11) 


~ 8 (~ 7) 


10 3 


314(159) 


3515 (2480) 


3312(2492) 


52 (27) 


311 (225) 


169(133) 


10 4 


381 (214) 


5963 (4488) 


7302 (5757) 


269 (123) 


2439(1618) 


1906 (1366) 


10 5 


368 (211) 


7258 (5575) 


10879 (8713) 


360 (197) 


4813 (3541) 


5191 (4023) 
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Late 


10 


< 1 « 1) 


~ 2 (~ 2) 


~ 1 (~ 1) 


< 1 « 1) 


< 1 « 1) 


< 1 « 1) 


10 2 


~ 5 (~ 4) 


42 (34) 


26 (22) 


< 1 « 1) 


~ 1 (~ 1) 


~ 1 (< 1) 


10 3 


195 (69) 


1132 (658) 


578 (400) 


~ 3 (~ 2) 


28 (23) 


18(15) 


10 4 


305 (152) 


3278 (2288) 


2987 (2229) 


172 (56) 


846 (493) 


390 (287) 


10 5 


380 (214) 


5752 (4314) 


6836 (5374) 


294 (143) 


3013 (2074) 


2629 (1939) 



Table 3: Number of FGST 5cr point sources in several WIMP annihilation scenarios for each of the 
five black hole masses and Early, Intermediate, and Late Zf, assuming both f DS = 1 and B f = l. 
In parentheses are the number of such sources with Galactic latitude |6| > 10°, to be compared 
with the 368 5a point sources with \b\ > 10° in the FGST First Source Catalog. 



the diffuse gamma ray flux from all faint spikes and make comparisons to the diffuse flux measured 
by FGST [HU. 

Since we use the diffuse flux to place limits on the population of black holes with DM spikes 
in the Milky Way halo, it is important that we not overestimate the diffuse flux. In other words, 
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we are careful to exclude from the analysis any spikes that could be identified as point sources in 
the one year FGST data. In Section |4~T1 we discuss the maximum distance D^ x for point sources 
detectable by FGST. Any black hole spikes located at distances less than D^ ax are excluded from 
our estimate of the diffuse flux. In fact, in the following analysis we take the more conservative 
approach of requiring an even larger minimal distance for spikes that contribute to the diffuse flux, 
Dmin', as discussed below. 

Ref. [57] estimate that the diffuse gamma-ray background produces 20 events per year per 
square degree above 1 GeV in FGST for \b\ > 60°, with a larger rate at lower Galactic latitudes. 
Here, we impose the condition that only spikes producing fewer than 20 events per year in FGST 
are included in the calculation of the diffuse gamma-ray flux. As mentioned above, we assume 
that a 5 a detection of a point source would require roughly 50 events. Since the number of events 
scales inversely with distance squared, we estimate that the minimal distance at which a DM spike 
would likely contribute to the diffuse flux is 



We note that sources located farther from our Solar System than Dj^ x and closer than D m ^' are 
included neither as point sources nor as part of the diffuse gamma-ray flux. 

In the following analysis, we assume that D^ ax and -D^fn ^ independent of Galactic lati- 
tude. However, in actuality, the background to detection of point sources does depend on Galactic 
latitude, and therefore both the flux required to generate a 5a excess and the flux below which the 
source would be too dim to be identified as a point source should also depend, to some degree, on 
the Galactic latitude. To be specific, at high Galactic latitudes the diffuse gamma-ray flux is lowest, 
so we expect that by choosing -D,^// as in Eq.[15l we may slightly overestimate the contribution of 
high |b| DM spikes to the diffuse gamma-ray flux. Sources resulting in 20 events per year above 1 
GeV in FGST for |6| > 60° represent a ~ 3ct fluctuation on the diffuse background, and therefore 
such sources may be bright enough to have been identified as point sources. However, Ref. ||59l 
find that 16(±1.8)% (with a systematic uncertainty of 10%) of the GeV isotropic diffuse back- 
ground is due to unresolved point sources, whose fluxes are large enough such that they could have 
been identified as point sources but have evaded detection due to selection effects. For sources at 
Galactic latitudes in the range 10° < |6| < 20°, 20 events per year above 1 GeV is only a ~ 2a 
fluctuation on the background. 

Given the assumption of a Galactic latitude-independent -D^f/' as explained above, we com- 
pute both the energy spectrum of the contribution to the diffuse gamma-ray background from all 
DM spikes at distances greater than -D^f/' and the angular distribution in the sky of the gamma-ray 
flux from those spikes. The latter quantity is for reference only, and may be sensitive to our choice 
of universal -D^f/ • Had we allowed -D^f/ t0 var Y w i m Galactic latitude, more spikes at low 
Galactic latitudes and fewer spikes at high Galactic latitudes would contribute to the diffuse flux, 
potentially resulting in a more anisotropic angular distribution than what has been calculated here. 
One should therefore view the anisotropics of the angular distributions presented here as minimal. 

In addition to neglecting the dependence on the Galactic latitude of the diffuse gamma-ray 
flux, we have also neglected any local deviations from the average value of the diffuse flux. In fact, 
the diffuse gamma-ray background is not entirely uniform, and local deviations from the average 




(15) 
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affect the brightness required of a point source for it to have been identified on top of the local 
background. We do not consider these issues here, and take instead a single value of D^f n for each 
model, applied to all regions of the sky. 

We see in the left panels of Figs. [6] and |7j that the all-sky averaged diffuse flux is well below the 
FGST measurement of the isotropic diffuse gamma-ray background in all models at low energies, 
but rises to the level of the FGST-measured flux at higher energies^. 

The diffuse flux from our models is expected to be relatively low in the following two cases: 
Case I: If the flux from individual DM spikes is very low such that even nearby spikes are not 
bright enough to be identified as point sources, we find that the total flux from all spikes is then 
also well below the measured diffuse flux. An example of a case which contributes little in the 
form of both point sources and the diffuse flux is the annihilation of heavy WIMPs in the DM 
spikes around very small black holes; Fig.HJillustrates the low luminosity of a single black hole in 
this example. 

Case II: Alternatively, if a significant fraction of the spikes in the Milky Way halo are bright enough 
to have been discovered as point sources, then only spikes very far from our Solar System would 
be faint enough to escape identification as point sources and therefore contribute to the diffuse 
gamma-ray flux. This is the case in many models where the individual point sources are very 
bright; when the black holes are very massive, when they formed very early, and if the dark matter 
is relatively light. Examples of this include annihilations to bb or W + W~ around black holes more 
massive than 1O 2 M , as shown in the left panels of Fig.[6l 

For scenarios that lie between the two extremes of Case I and Case II, the expected diffuse 
flux of gamma-rays from all unresolved DM spikes in the Milky Way halo may be significantly 
larger than the FGST-measured diffuse flux. Therefore, for each black hole mass, Zf, and dark 
matter mass and annihilation mode, a constraint on f DS may be derived from the requirement that 
the diffuse flux of gamma-rays from DM spikes not produce a significant excess over the FGST 
observation. These constraints will be discussed in Section[5] 

It is also interesting to examine the anisotropy of the expected diffuse gamma-ray flux, namely 
to study the number of events as a function of angle with respect to the Solar System- Galactic 
center axis. This anisotropy is plotted in the right panels of Figs. [6] and [7J Previously it was sug- 
gested that the anisotropy might be used to discriminate among models of Milky Way substructure 
and dark matter annihilations Il60l . Though such an analysis has not yet been done with data from 
FGST, here we consider the behavior of the anisotropy as an indication of whether DM spikes near 
the Galactic center contribute to the diffuse gamma-ray flux. 

For very low luminosity spikes, as in Case I, few, if any, spikes are bright enough to have been 
identified as point sources. If this is the case, then nearly all spikes contribute to the diffuse signal, 
including those near the Galactic center. The diffuse gamma-ray flux would therefore be peaked 
towards to Galactic center. As one can see in the right panels of Figs. [6] and [71 for most dark 
matter models shown, the spikes surrounding 1O 2 M black holes (thick curves) are indeed so dim 
that those from the Galactic center region contribute to the diffuse flux. The exceptions to this are 
the two most luminous dark matter models shown, 6100 and W 100, which exhibit more isotropic 

4 For the sake of the clarity of Figs. [6] and [TJ we have plotted only two choices of central black hole mass, 10 2 and 
10 3 Af Q . Larger central black holes may result in either larger or smaller gamma-ray fluxes, depending on the WIMP 
annihilation model, generally not differing by more than an order of magnitude from the uibh = 1O 2 M case. 
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E (GeV) cos(0) 

Figure 6: Diffuse gamma-ray flux from DM spikes for models blOO, blT, W100 and WIT, from 
top to bottom. In each panel the Early (green, dotted), Intermediate (red, solid), and Late (blue, 
dashed) formation scenarios are shown for central black holes of mass 1O 2 M (thick curves) and 
1O 3 M (thin curves). In the left panels, the all-sky average flux is shown, with the FGST-measured 
EGB as black points 129]. In the right panels, the maximal anisotropy is displayed. 
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Figure 7: Diffuse gamma-ray flux from DM spikes for models rlOO, rlT, jilOO and filT, from top 
to bottom. Curves are as described in Fig\$ 



signals. 

All other factors being equal, the DM spikes surrounding 10 3 M Q black holes are nearly an 
order of magnitude brighter than those surrounding 10 2 M Q black holes. These higher-luminosity 



19 



DM spikes are generally bright enough that the closest ones must be at least as far from our Solar 
System as the Galactic center, as evidenced by the lack of peaking toward the Galactic center for 
the thin curves in the right panels of Figs. [6] and [71 The exceptions that exhibit a peak towards 
the Galactic center for 1O 3 M black holes are models /ilT and rlT, the least luminous models 
considered here. For even larger black holes, not shown in Figs. [6] and [/J the luminosity of each 
spike may be so high that only spikes in the very outer regions of the Milky Way halo (if any) 
contribute to the diffuse gamma-ray flux, as in Case II. This results in a nearly-completely isotropic 
diffuse flux. 

In contrast to the magnitude of the diffuse flux, which would be small in both Case I and Case 
II, the anisotropy is expected to be different in the two cases: Case I would result in some peaking 
of the signal towards the Galactic center, while the signal from Case II would be nearly isotropic. 
However, it must be pointed out that any anisotropy would be most evident near the Galactic 
center, where there are many baryonic gamma-ray sources, making the diffuse component difficult 
to disentangle, and where the survival of DM spikes is less certain. If spikes in the central region 
were disrupted, any anisotropy may have been washed out. Therefore, if there is a peak in the 
angular distribution of the diffuse gamma-ray flux towards the Galactic center, it may be evidence 
of a Case /-like population of DM spikes, but the lack of a strong anisotropy is not necessarily 
evidence for a Case //-like population. 

In all scenarios discussed above, the star formation history plays a significant role in so far as 
the total number of DM spikes in the Milky Way halo is greatest if star formation persisted to low 
redshift. If Population III. 1 star formation ended at high redshift, as for Early z/, there would be 
significantly fewer black hole remnants today in the Milky Way halo. This scenario, identified by 
the green curves in Figs. |6] and |7J universally results in a lower diffuse flux than the other scenarios 
examined here. The Intermediate and Late Zf scenarios, identified by the red and blue curves, 
respectively, typically have results that are quite similar to each other, both for the total gamma-ray 
fluxes and for the spectra. Again, we emphasize that we have taken f DS = 1 in these calculations. 

5 Constraining fns 

To this point we have focused mainly on f DS = 1, though it is possible that f DS <C 1. It is 
possible to constrain f DS in two ways: First, from the FGST measurement of the diffuse gamma- 
ray background, we require that the diffuse flux from dark matter annihilations around spikes in 
the Milky Way halo not exceed the measured flux in any of the nine FGST energy bins in Ref. [|29l 
by more than 3a. The diffuse gamma-ray flux from annihilations to final state / may be rewritten 
for f DS ^ 1 as 

$/M = /dsX%(/ ds = 1). (16) 

One can then extract an upper limit on f DS from the diffuse background for each choice of Zf and 
central black hole mass in each of the dark matter annihilation models. The resulting maximal 
values of fos are presented in Fig. [8] for each of the WIMP annihilation models in Table |2] as the 
open points, with Early, Intermediate, and Late Zf represented in green, red, and blue, respectively. 
Note that some scenarios are entirely unconstrained (maximal f DS = 1), in which case the markers 
are superposed on each other. These constraints are strongest for scenarios in which the intrinsic 
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luminosity of individual DM spikes is moderate; that is, bright enough that the signal is strong, but 
not so bright that many /all spikes in the Milky Way halo are seen as point sources. The constraints 
are also weakest in all cases for Early termination of star formation, as there are fewer potential 
nearby DM spikes today. 

The second way to constrain fns is by requiring that one should expect to find less than one 
DM spike within D^f n , the minimal distance at which a spike may be located such that it is not 
brighter than the brightest point source in the FGST First Source Catalog. In this case we rewrite 
the number density of DM spikes as a function of distance from the Galactic center, N sp (R), as 



The maximum f^s as found in this method is displayed in Fig. [8] as the solid points in each 
scenario. These are the most conservative limits on fas, as they are based on the flux not exceeding 
that of the brightest FGST point source, the Vela pulsar. If we require that the flux from a DM 
spike not exceed that of the brightest unassociated point source, stronger limits on f^s would be 
obtained, though there is no gaurantee that there is not a DM spike along our line of sight to Vela. 

We see from the solid points in Fig. [8] that the point source constraint on f DS is strongest in 
all cases for large central black hole masses and low WIMP masses, when the intrinsic luminosity 
of each individual spike is largest. Again, they are weak in the case of Early termination of star 
formation, as there are fewer black holes expected in the Milky Way halo, and therefore fewer in 
the vicinity of our Solar System. 

To this point we have not properly addressed the uncertainty in the DM spike distribution near 
the Galactic center. In fact, in the central region of the Milky Way halo, subhalo evolution is 
affected by processes such as dynamical friction and tidal mass loss lloTll . resulting in a dearth of 
dark matter substructures at small Galactic radii. However, as pointed out in Ref. [36J, black holes 
and their surrounding dark matter spikes can survive tidal disruption even near the Galactic center. 
Still, Ref. 11361 finds that the distribution of DM spikes is very uncertain, and potentially even zero, 
in the inner ~ 3 kpc (see, for example, Fig. 1 of (33). 

Given this uncertainty, we calculate also the constraints on fos m each scenario, neglecting all 
DM spikes in the inner 5 kpc of the Galaxy. In Fig.|9l these results are presented. If fos = 1, for 
Early, Intermediate, and Late Zf, we expect roughly 154, 407, and 895 DM spikes in the inner 5 
kpc, corresponding to 38%, 10%, and 4% of the total number of DM spikes in the Milky Way halo, 
respectively. The constraints are strikingly similar to those in Fig. [8] In general, the limits on fos 
from point source brightness come from much smaller distances, as can be seen in Fig.[5l while the 
limits from the diffuse flux are sensitive primarily to the bulk of the halo at large Galactic radii. As 
a result, there are only a few cases in which the constraints visibly differ, all of which concern the 
nearest point source in scenarios where individual DM spikes are bright enough that D^f n > 8.5 
kpc. With the exception of these few cases, our constraints on f DS are robust with respect to the 
distribution of DM spikes near the Galactic center. 

Finally, we return to the issue of the maximum possible brightness of DM spikes. In the pre- 
ceeding analysis, we consider only the most conservative criterion for the maximum DM spike 



N SP (R, Ids) = fos x N sp (R, f DS = 1), 
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and find f DS such that 
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Figure 8: Maximum fps as a function of central black hole mass for the WIMP annihilation models 
in Table\2\ Green circles, red squares, and blue diamonds are for Early, Intermediate, and Late Zf, 
respectively. The solid markers are the limits from point source brightness, while the open markers 
are from the diffuse gamma-ray flux, as described in the text. In each panel, the branching fraction 
to the relevant final state Bf = 1. Note that the range of fas displayed differs from panel to panel. 
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Figure 9: Maximum f^g as in Fig. \8\ but excluding the inner 5 kpcfrom the Galactic center. 

brightness, namely that the brightest DM spike must not be brighter than the brightest FGST point 
source, Vela. Given the size of Vela in the sky, it is highly unlikely that the brightest DM spike 
happens to lie along our line of sight to Vela. In fact, even if it does, it is then extremely unlikely 
that the second brightest DM spike would also lie along our line of sight to another bright asso- 
ciated FGST source (and so on for subsequently dimmer DM spikes). It is much more likely that 
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Figure 10: Maximum fjjg as in Figs. \t^and® excluding the inner 5 kpc from the Galactic cen- 
ter and with the maximum fos from point source brightness determined by the brightness of the 
brightest unassociated point source. The maximum f^s from the diffuse gamma-ray flux is as in 
Fig. [9] The range of fos displayed is different than in previous figures for p and r final states. 

the brightest DM spike does not lie precisely along our line of sight to Vela, or any other of the 
brightest FGST point sources. 
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Given the unlikelihood that DM spikes are hiding along our line of sight to the brightest asso- 
ciated FGST point sources, here we explore the possibility that the brightest DM spike must not be 
brighter than the brightest unassociated point source, which has an energy-integrated gamma-ray 
flux of 5.78 x 1CT 8 cm^s -1 , approximately 1 /22 that of Vela. In Fig. [TOl we present the maxi- 
mum value of fos as determined by Eq. [18] with the more stringent maximal flux. As in Fig. 
we exclude any spikes located within 5 kpc of the Galactic center. Note that the ranges of fns 
displayed for ji and r final states are different from those in Figs. [8] and [9l and that there has been 
no change to the constraint derived from the diffuse gamma-ray flux. It is clear that the nearest 
point source must now be even farther away, leading to much stronger constraints on fos- Again, 
the limits on fns are strongest for the most luminous individual spike scenarios (low WIMP mass, 
annihilations to bb or W + W~), and for the star formation histories that result in the largest number 
of DM spikes near our Solar System. For annihilations of 100 GeV WIMPs to bb or W + W~, there 
are significant limits on f DS even for 1OOM black holes, as would be expected from standard 
Population III. 1 stars in the absence of a DS phase. 

The limits on f DS presented here have all been obtained assuming, as discussed in Sec. 14. 11 that 
the flux from an individual spike necessary for it to have been identified as an FGST point source 
does not depend on the location of the spike in the sky. In fact, the flux required for a point source 
to have been identified does depend on the diffuse gamma-ray emission near the spike. However, 
since some spikes would have been easier to identify while others may be in regions of the sky that 
inhibit identification, we estimate that our results are not sensitive to deviations from the average 
flux we've used. Specifically, since D^ n only depends on the flux from the brightest point source, 
the constraint from near point sources is not affected at all, while the constraint from the diffuse 
flux comes from an ensemble of many sources, and will therefore be minimally affected. 

If, however, the number of events required to identify a point source is significantly larger or 
smaller than the 50 per year that we have assumed, our results would change somewhat. In either 
case, there is no difference in the constraint arising from the brightest spike. If the number of events 
required is larger than 50 per year, more spikes would not be bright enough to have been identified 
as point sources, so the diffuse gamma-ray flux from all faint (non-point source) spikes would be 
larger than that presented here, and therefore the constraint on f DS arising from the diffuse flux 
would be stronger than what is presented in Figs. |8lfT0l Conversely, if sources producing fewer 
than 50 events per year in FGST have been identified as point sources, then fewer spikes would 
contribute to the diffuse flux, and the constraint on f DS arising from the diffuse flux would be 
weaker than what is presented in Figs. l8~lfT0l Again, for the calculation of the diffuse gamma-ray 
flux we include only spikes that result in less than 20 events per year above 1 GeV. If more (fewer) 
spikes actually contribute to the diffuse flux, then our constraint on f DS from the diffuse flux would 
improve (weaken), while our constraint from nearby point sources would be unaffected^]. 

5 We would also expect that there would be fewer (more) spikes that should have appeared as point sources in the 
First Source Catalog, and therefore fewer (more) spikes in Table[3] 



25 



6 Antimatter and Neutrino Signatures of Nearby DM Spikes 



In general, dark matter annihilations result in other potentially-detectable products in addition to 
photons. For example, neutrinos, electrons and positrons, and a host of other Standard Model final 
states may be produced, depending on the annihilation mode. Neutrinos and antimatter present 
very interesting and potentially fruitful avenues in which to search for dark matter annihilations in 
nearby structures. While there are currently only limits on the dark matter annihilation rate to neu- 
trinos and antiprotons, several recent experiments sensitive to cosmic ray electrons and positrons 
have found anomalous signals. 

6.1 Positrons: Comparison with PAMELA data 

The Payload for Antimatter Matter Exploration and Light-Nuclei Astrophysics (PAMELA) has 
observed an excess of cosmic ray positrons (relative to electrons) between 10 and 100 GeV ll62l . 
confirming the hints from earlier experiments such as the High Energy Antimatter Telescope 
(HEAT) [63 J and the AntiMatter Spectrometer (AMS-01) ll64l . A surplus of cosmic ray electrons 
and/or positrons has also been confirmed by FGST [65 J, however the spectral feature is less pro- 
nounced than that previously reported by the Advanced Thin Ionization Calorimeter (ATIC) ll66l . 

Although an excess of cosmic ray positrons may come from more conventional astrophysical 
sources such as pulsars 11671 , there has been a great deal of interest in the possibility that these 
signals might be a consequence of dark matter particles annihilating in the local halo of the Milky 
Way. Efforts to produce such signals with dark matter, however, have faced two major challenges. 
First, a very large annihilation rate is required, much larger than that expected from a smoothly- 
distributed thermal relic annihilating in the Milky Way halo. Second, given the large annihilation 
rate, for most annihilation modes a positron excess would be accompanied by an excess of antipro- 
tons. However no such excess of cosmic ray antiprotons has been observed. 

These two challenges have been addressed in a variety of ways, including string theory moti- 
vated supersymmetric models with non-thermal wino dark matter, as proposed by 11681 . Of partic- 
ular relevance to our current work is Ref. 11691 , which suggests that annihilations in a nearby clump 
of dark matter may be responsible for the positron excess. It may be extremely unlikely that a 
DM subhalo with the required luminosity and distance exists H48U70II , but the required luminosities 
and distances do lie within the range of DM spike scenarios considered here and in |[36] - [39l . For 
example, an 800 GeV WIMP annihilating to W + W~ in the DM spike surrounding a 1O 2 M black 
hole has a luminosity of ~ 10 37 7/s, according to Fig. @] It was found by Ref. |[69l that this is 
roughly the luminosity required to explain the PAMELA positron excess if the dark matter clump 
or spike is located approximately 1 kpc from our Solar System. From the left panel of our Fig.[5J 
we see that for a 10 2 M Q black hole and annihilations to W + W~, the nearest spike as constrained 
by the FGST point source data may be as close as ~ 1/2 kpc, potentially over-producing cosmic 
ray positrons. If that DM spike is 1 kpc from our Solar System, in agreement with PAMELA, it 
may already appear in the FGST catalog. 

The feasibility of explaining the excess of cosmic ray positrons with nearby DM spikes sur- 
rounding 1O 5 M black holes was also addressed in Ref. [[391 , where it was confirmed that a DM 
spike ~ 1 kpc from our Solar System could indeed reproduce the PAMELA excess. Moreover, it 
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was shown that a range of annihilation rates and spike distances could explain the PAMELA excess 
without overproducing the combination of electrons + positrons observed by FGST. The implica- 
tions of the existence of more than one nearby DM spike may be important as measurements of 
the flux of cosmic ray electrons and positrons improve. A more detailed analysis of the interplay 
between gamma-ray and antimatter constraints in these scenarios, including the effects of multiple 
spikes according to the local distribution as predicted by VL-II, we leave for future work 11711 . 

6.2 Neutrinos 

Although the photon flux generally imposes a stronger constraint on dark matter annihilations, the 
neutrino flux may become more significant if there is a substantial branching fraction to leptonic 
final states. Searches for a diffuse neutrino flux from dark matter annihilations are limited by the 
appreciable background of atmospheric neutrinos above ~ 40 MeV, though it has been shown that 
if there is a substantial branching fraction to leptonic final states, the neutrino flux may be a very 
useful tool to learn about the properties of dark matter through its annihilations in the Milky Way 
halo 11721 - 1741 and in the dark matter halos of Milky Way dwarf satellite galaxies [17511761 . 

The prospects for detecting neutrinos from dark matter annihilations in the spikes surrounding 
1O 5 M black holes were investigated in Ref. [77], where the number of DM spikes producing suf- 
ficient event rates at several neutrino telescopes was calculated. Ref. ir77T| concludes that, based on 
the distribution of such black holes from [36J, neutrinos from 1O 5 M black holes could be detected 
by both the IceCube neutrino detector at the South Pole and ANTARES, in the Mediterranean Sea. 

The most stringent limits on the flux from point sources in the Northern Hemisphere have 
been set by the AMANDA-II Collaboration [78], though these limits apply only to neutrinos with 
E v > 1.6 TeV, a threshold too large to be relevant for the models considered here. The full 
IceCube/DeepCore neutrino detector will have improved sensitivity for point source searches with 
neutrino energies potentially as low as ~ 10 GeV [1791 . The ANTARES neutrino telescope claims 
the best limits on the neutrino flux from point sources in the Southern Hemisphere for 10 GeV 
<E U < 100 TeV of Eld$/dE v < few x 1(T 7 GeVcm~ 2 s^ 1 |[80l. 

Super- Kamiokande (SK), located in the Mozumi Mine in Japan, is sensitive to neutrinos from 
the Southern Hemisphere down to E v = 4.5 MeV, where the flux of solar neutrinos dominates. 
Recently, SK performed a search for neutrino point sources and found that $ < 10~ 7 cm _2 s _1 for 
E v > 1.6 GeV at declinations between -90° and ~ 50° MM- We plot in Fig. \TT\ the minimum 
distance at which a DM spike must be located in order not to exceed this flux, assuming the source 
is located in the field of view of SK, as in Fig. |5l Here we use full three flavor vacuum oscilla- 
tions [82], and consider only neutrinos that arrive at Earth as u^, as in Ref. lUffll . Comparing Fig.fTTl 
and the left panel of Fig. [51 we see that neutrino point source brightness provides a slightly stronger 
constraint on the minimal distance of the nearest spike in all scenarios shown if the brightest DM 
spike must not be brighter than Vela in gamma-rays. 

Of course, unjustified assumptions have been made here; most obviously that the SK limit on 
the flux from point sources applies to the full sky. There are many caveats, as well, such as the fact 
that the limit was placed under the assumption that all point sources emit neutrinos with a spectral 
index of 7 = 2, such that oc <E>^ ' E^ 1 '. For softer spectra (larger 7), the limits are considerably 
weaker [1811 . For these reasons we abandon, for now, the possibility of using neutrinos to constrain 
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Figure 1 1 : The minimum distance at which a spike may be found such that the flux of neutrinos 
from dark matter annihilation in the spike does not exceed the brightest neutrino point source 
as measured by Super Kamiokande. This distance is plotted as a function of WIMP mass for 
Intermediate Zf. From top to bottom, the contours are for models blOO (blue dashed), W100 
(orange dashed), rlOO (magenta dashed), jilOO (green dashed), WIT (orange solid), b IT (blue 
solid), rlT (magenta), and filT (green solid). The black horizontal line indicates our distance 
from the GC. 

dark matter annihilations in density spikes surrounding 10 — 1O 5 M black holes. We note, however, 
that future neutrino observations with the DeepCore supplement to the IceCube detector, which 
will probe energies down to E v ~ 10 GeV (possibly with full-sky coverage) [|83l . as well as 
future observations from SK, may provide possibilities for novel signatures of nearby DM spikes 
in neutrinos. The interplay of neutrino and gamma-ray astronomy in the indirect detection and 
identification of dark matter may be very important in the near future. 

7 Discussion and Conclusions 

In this paper we have examined 10 — 1O 5 M black holes with DM spikes that formed in early mini- 
halos and still exist in our Milky Way Galaxy today. Such black holes may arise as the remnants 
of Dark Stars after the dark matter fuel is exhausted and thermonuclear burning runs its course. 
Since the redshift at which the increasing UV background and/or metal enrichment results in the 
truncation of Pop. III. 1 star formation is poorly constrained, we have examined three scenarios: 
Early (zf ~ 23), Intermediate (zj ~ 15), and Late (zf ss 11). We determined the z = distri- 
bution of DM spikes throughout the Galactic halo from the Via Lactea II cosmological N-body 
simulation [|48]|. We find 409, 7983, and 12416 DM spikes in the Milky Way for the Early, Inter- 
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mediate, and Late scenarios, respectively, for the case where fns = 1 (every minihalo hosts a BH 
with a spike, no mergers). As the DM spikes around the black holes are the sites of significant DM 
annihilation, we examined the signatures of the spikes in gamma-rays, and commented on searches 
in e + /e _ and neutrinos. 

The focus of this paper is the gamma-rays from DM annihilation in the spikes around black 
holes in the Galaxy, and comparison with data from FGST. We find that some or all of the unas- 
sociated point sources observed by FGST could be black holes with DM spikes. Indeed, Table 2 
shows the estimated number of point sources in the Milky Way halo for a variety of DM models, 
and one can see that in most cases more than the 368 observed sources are present. It is exciting to 
imagine that some significant fraction of the point sources might be due to DM annihilation near 
black holes in our Galaxy. For further discussion of DM point sources in the FGST First Source 
Catalog, see Ref. 1571 . 

Additionally, one can use the FGST observations to place limits on the properties of DM spikes 
around black holes. We compare the gamma-ray flux from our models to both that from point 
sources and the diffuse flux observed by FGST, and ensure that they do not produce gamma-rays 
in excess of what is observed. We have found the maximum fraction of Pop. III. 1 star-forming 
minihalos that could have contained DSs and their black hole remnants. Our results are shown in 
Figs. [8lfT0l In general, it is clear that the bounds are the strongest for the largest black hole masses 
and if star formation persisted to low redshift (i.e. the Late model where star formation ends at 
zj = 11). We have also found that our bounds are robust with respect to uncertainties related to 
dynamics near the Galactic center, in that neglecting all DM spikes in the inner 5 kpc of the Galaxy 
does not substantially affect the limits. 

We note that our limits depend very sensitively on the DM annihilation channel and on the 
black hole mass. One can see that even the conservative constraints from bright point sources, 
shown in Figs. [8] and [9l are stronger than those from the diffuse flux in the case of annihilation 
primarily to bb or W + W~ . These cases are most typical of Minimal Standard Supersymmetric 
Model (MSSM) neutralinos. The bounds become more restrictive for higher mass black holes. For 
example, the upper left panels in both Figs. [8] and [9] illustrate the bounds for 100 GeV WIMPs 
with annihilation primarily to bb. There is essentially no constraint for m BH < 100M Q , while the 
constraint on the fraction of minihalos containing black hole spikes becomes f DS B b i < 1/10 for 
IOOOMq black holes for Intermediate and Late Zf. For the case of 100 GeV WIMPs and even more 
massive 1O 5 M black holes, the bound becomes fosB^i < .01, regardless of Zf. If the number of 
DSs with these properties is low, it becomes more difficult to find them in upcoming observations 
with James Webb Space Telescope (JWST). On the other hand, for TeV-mass WIMPs, the upper 
right hand panel of Fig. [9] shows a weaker constraint for the 1O 5 M black holes of f DS B b i < 1/10. 
Thus one out of every ten minihalos with the right properties for Pop. III. 1 star formation might 
have hosted a DS with these characteristics. This is a good case for discovery of DS in JWST J9]|. 
From the remaining panels in Figs. [8] and [9l one can see that the constraints are relatively weak for 
WIMPs that annihilate to leptonic final states. 

While scenarios involving the larger black holes considered here are decidedly constrained for 
all models in Figs.[8]and|9l it is remarkable that we are able to find conservative limits on f DS even 
for IOOMq black holes for some dark matter models. Black holes of this size are roughly what is 
expected from standard Pop. III. 1 star formation, if there is no DS phase. We also point out that 
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the smallest black holes we consider here, with tubh = 10M Q , are actually disfavored by both 
theory ll84l and simulations J2l[85l as the remnants of Pop. III. 1 stars, which are expected to be 
considerably larger. 

In addition to considering the most conservative case that allows for the possibility that the 
brightest DM spike could be located along our line of sight to the brightest FGST point source, we 
also consider the more likely scenario that the brightest DM spike is in fact not hiding along our line 
of sight to any of the brightest sources in the FGST catalog that are associated with astrophysical 
sources in other wavelengths. If the brightest DM spike must not be brighter than the brightest 
unassociated object in the FGST First Source Catalog, fog must indeed be quite small for tubh > 
1OOOM for most dark matter models explored. For annihilations of 100 GeV WIMPs to bb or 
W + W~, the limits on f DS are quite significant even for 100Af o black holes, as would be expected 
from standard Population III. 1 stars in the complete absence of a DS phase. In order to trust these 
stronger limits, however, we must be completely confident in our understanding of the brightest 
associated FGST point sources and confident that there is no DM spike hiding along the line of 
sight to any of them. 

We stress that the limits derived here on fog are related to the initial fraction of minihalos 
in which Population III. 1 stars formed and the fraction of DM spikes subsequently disrupted by 
mergers as fos = /z?s(l ~ f merged)- While our estimates show that mergers don't change the 
number of DM spikes in our Galaxy by more than a factor of 2, a detailed understanding of black 
hole mergers is required to translate our constraints on fog to constraints on the fraction of mini- 
halos that were capable of hosting Population III. 1 stars that actually did, f^ s . This information 
could then be used to place an explicit limit on the Population III. 1 star formation rate. If/when 
the annihilation properties of particle dark matter become known, we may be provided with some 
interesting hints about the formation of the first stars. 

Here we have restricted our work to black holes with m BH < 1O 5 M and WIMP masses in 
the range 100 GeV - 2 TeV. Larger DSs are indeed possible [j9j|, the most massive of which are 
expected to be sufficiently bright to be detected by the Hubble Space Telescope (HST) ll24l . In 
fact, Ref. Il24l uses HST data to constrain the formation rate and lifetime of 1O 7 M SMDSs. It 
would be interesting to directly compare the constraints derived from FGST data with those from 
HST. It would also be very interesting to examine the possibility of lighter O(10) GeV WIMPs. On 
the one hand, the annihilation rate (which scales inversely with the WIMP mass) would be larger, 
so the sources would be brighter. On the other hand, the photon energies would be lower so it's 
possible that these light WIMP scenarios are already tightly constrainted by the EGRET-measured 
photon flux. This would be an interesting case to study, particularly since current data from direct 
detection experiments seem to favor this WIMP mass range; DAMA/LIBRA and CoGeNT data 
may be consistent with such a value |[86] - [891 , though recent reports may indicate incompatibility 
with XENON 10 data 1901 . 

We also briefly discussed the relevance of DM spikes in relation to the anomalous cosmic ray 
positron excess measured by PAMELA. DM spikes seem to be consistent with the dark matter 
annihilation interpretation of the excess. Further study of the distribution of spikes in the Milky 
Way halo and compatibility with cosmic ray antimatter spectra is merited [1711 . 

In this paper, we have shown how dark matter structures that exist today carry information 
about the properties of the first stars in the Universe, and we have explored a few ways to learn 
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about them through existing gamma-ray data from FGST. We look forward to the day when the 
properties of particle dark matter reveal themselves, and the implications for Dark Stars and dark 
matter astronomy unfold. 

Acknowledgements 

J.D. is supported by the Swiss National Science Foundation. K.F. thanks the Department of Energy 
and the Michigan Center for Theoretical Physics for support, and the Aspen Center for Physics for 
hospitality during the course of this research. K.F. thanks Paul Shapiro for helpful conversations. 
PS. is supported by the National Science Foundation under Grant No. PHY-0455649, and thanks 
the Michigan Center for Theoretical Physics for providing a stimulating working environment 
during the Dark Stars Workshop. PS. also thanks Eiichiro Komatsu for helpful conversations. 
D.S. is supported by the Department of Energy. 

References 

[1] Z. Haiman, A. A. Thoul and A. Loeb, Astrophys. J. 464, 523 (1996) 
|arXiv:astro-ph/9507TTTl . 

[2] V. Bromm, P. S. Coppi and R. B. Larson, Astrophys. J. 527, L5 (1999) 
[ |arXiv:astro-ph/9910224| . 

[3] V. Bromm and R. B. Larson, Ann. Rev. Astron. Astrophys. 42, 79 (2004) 
[ |arXiv:astro-ph/031 10T9| . 

[4] R. Barkana and A. Loeb, Phys. Rept. 349, 125 (2001) |arXiv : astro-ph/00 10468Q . 

[5] E. Ripamonti and T. Abel, ||arXiv: astro-ph/050 7T30[ 

[6] N. Yoshida, T. Abel, L. Hernquist and N. Sugiyama, Astrophys. J. 592, 645 (2003) 
narXiv:astro-ph/0301645l . 

[7] D. Spolyar, K. Freese and P. Gondolo, Phys. Rev. Lett. 100, 05 1 101 (2008) HarXiv:0705.0521l 
[astro-ph]]. 

[8] K. Freese, P. Bodenheimer, D. Spolyar and P. Gondolo, Astrophys. J. 685, L101 (2008) 
HarXiv:0806.0617l [astro-ph]] . 

[9] K. Freese, C. Hie, D. Spolyar, M. Valluri and P. Bodenheimer, larXiv: 1002.22331 [astro- 
ph.CO]. 

[10] F. Iocco, Astrophys. J. 677, LI (2008) HarXiv:0802.0941i rastro-phll. 

[11] K. Freese, D. Spolyar and A. Aguirre, JCAP 0811, 014 (2008) I arXiv:0802. 17241 rastro-phll . 



31 



[12] K. Freese, P. Gondolo, J. A. Sellwood and D. Spolyar, Astrophys. J. 693, 1563 (2009) 
HarXiv:0805.3540l [astro-ph]] . 

[13] M. Taoso, G. Bertone, G. Meynet and S. Ekstrom, PoS IDM2008, 076 (2008). 

[14] S. C. Yoon, F. Iocco and S. Akiyama, Astrophys. J. 688, LI (2008) HarXiv:0806.2662l fastro- 
ph]]. 

[15] E. Ripamonti, F. Iocco, A. Bressan, R. Schneider, A. Ferrara and P. Marigo, PoS IDM2008, 
075 (2009) HarXiv:0903.0346l rastro-ph.CQll . 

[16] F. Iocco, A. Bressan, E. Ripamonti, R. Schneider, A. Ferrara and P. Marigo, Mon. Not. Roy. 
Astron. Soc. 390, 1655 (2008) HarXiv:0805.4016l fastro-phll . 

[17] D. Spolyar, P. Bodenheimer, K. Freese and P. Gondolo, Astrophys. J. 705, 1031 (2009) 
HarXiv:0903.3070i rastro-ph.CQll. 

[18] D. R. G. Schleicher, R. Banerjee and R. S. Klesser, Phys. Rev. D 78, 083005 (2008) 
HarXiv:0807.3802l [astro-ph] ] . 

[19] D. R. G. Schleicher, R. Banerjee and R. S. Klessen, Phys. Rev. D 79, 043510 (2009) 
HarXiv:0809.1519l [astro-ph]] . 

[20] E. Ripamonti, F. Iocco, A. Ferrara, R. Schneider, A. Bressan and P. Marigo, arXiv: 1003.0676 
[astro-ph.CO]. 

[21] H. Umeda, N. Yoshida, K. Nomoto, S. Tsuruta, M. Sasaki and T. Ohkubo, JCAP 0908, 024 
(2009) HarXiv:0908.0573l [astro-ph.COll. 

[22] S. Sivertsson and P. Gondolo. larXiv: 1006.00251 rastro-ph.CQl . 

[23] E. Zackrisson etal, Astrophys. J. 717, 257 (2010) llarXiv: 1002.33681 rastro-ph.CQll . 

[24] E. Zackrisson et al. , larXiv: 1006.048T1 [astro-ph.CO] . 

[25] J. C. Tan, IAU Conf. Proc. 255, 24 (2008) larXiv:0808.3918l lastro-phl. 

[26] A. Heger and S. E. Woosley, Astrophys. J. 567, 532 (2002) |arXiv:astro-ph/0107037| . 

[27] A. Heger, C. L. Fryer, S. E. Woosley, N. Langer and D. H. Hartmann, Astrophys. J. 591, 288 
(2003) [arXiv:astro-ph/0212469| . 

[28] T. L. Collaboration. larXiv: 1002.22801 lastro-ph.HEl . 

[29] A. A. Abdo et al. [The Fermi-LAT collaboration], Phys. Rev. Lett. 104, 101101 (2010) 
llarXiv: 1002.36031 rastro-ph.HEll . 

[30] P. Gondolo and J. Silk, Phys. Rev. Lett. 83, 1719 (1999) [ 'arXiv:astro-ph/990639 1| . 



32 



[31] P. Ullio, H. Zhao and M. Kamionkowski, Phys. Rev. D 64, 043504 (2001) 
garXiv:astro-ph/010f48Tl . 

[32] M. Preto, D. Merritt and R. Spurzem, Astrophys. J. 613, L109 (2004) 
garXiv:astro-ph/0406324) . 

[33] G. Bertone and D. Merritt, Phys. Rev. D 72, 103502 (2005) [arXiv:astro-ph/0501555]. 



[34] G. Bertone and D. Merritt, Mod. Phys. Lett. A 20, 1021 (2005) [arXiv:astro-ph/0504422] 



[35] H. S. Zhao and J. Silk, Phys. Rev. Lett. 95, 01 1301 (2005) | [arXiv:as tro-ph/0 5016251 1. 

[36] G. Bertone, A. R. Zentner and J. Silk, Phys. Rev. D 72, 103517 (2005) 
garXiv:astro-ph/0509365| . 

[37] G. Bertone, M. Fornasa, M. Taoso and A. R. Zentner, New J. Phys. 11, 105016 (2009) 
HarXiv:0905.4736l [astro-ph.HE]] . 

[38] M. Taoso, S. Ando, G. Bertone and S. Profumo, Phys. Rev. D 79, 043521 (2009) 
HarXiv:08 11.44931 [astro-ph]] . 

[39] T. Bringmann, J. Lavalle and P. Salati, Phys. Rev. Lett. 103, 161301 (2009) HarXiv:0902.3665l 
[astro-ph.CO]]. 



[40] Z. Haiman, T. Abel and M. J. Rees, Astrophys. J. 534, 11 (2000) [arXiv:astro-ph/9903336]. 



[41] K. Ahn and P. R. Shapiro, Mon. Not. Roy. Astron. Soc. 375, 881 (2007) [AIP Conf. Proc. 
990, 377 (2008)] [ |arXiv:astro-ph/0607642[ l. 

[42] T. H. Greif and V. Bromm, Mon. Not. Roy. Astron. Soc. 373, 128 (2006) 
HarXiv:astro-ph/0604367| . 

[43] J. H. Wise and T. Abel, Astrophys. J. 685, 40 (2008) larXiv:0710.3 1601 fastro-phl . 

[44] C. F. McKee and J. C. Tan, Astrophys. J. 681, 771 (2008) larXiv:071 1.13771 fastro-phl. 

[45] M. Sakuma and H. Susa, Astrophys. J. 698, 155 (2009) HarXiv:0904.2355i rastro-ph.COll. 

[46] M. Trend and M. Stiavelli, Astrophys. J. 694, 879 (2009) [ arXi v:090 1 .07TT1 rastro-ph.COll . 



[47] J. H. Wise and T. Abel, Astrophys. J. 684, 1 (2008) arXiv:07 10 .4328 [astro-ph]. 

[48] J. Diemand, M. Kuhlen, P. Madau, M. Zemp, B. Moore, D. Potter and J. Stadel, Nature 454, 
735 (2008) HarXiv:0805. 12441 rastro-phn. 

[49] A. A. Dutton, F. C. van den Bosch, A. Dekel and S. Courteau, Astrophys. J. 654, 27 (2006) 
garXiv:astro-ph/0604553l . 

[50] J. Binney and S. Tremaine, Galactic Dynamics: Second Edition, Princeton University Press, 
Princeton, New Jersey (2008). 

33 



[51] J. F. Navarro, C. S. Frenk and S. D. M. White, Astrophys. J. 462, 563 (1996) 
garXiv:astro-ph7950 80251 . 

[52] C. Hie, K. Freese and D. Spolvar. larXiv: 1008.03481 rastro-ph.COl . 

[53] G. R. Blumenthal, S. M. Faber, R. Flores and J. R. Primack, Astrophys. J. 301, 27 (1986). 

[54] D. Merritt and M. Valluri, Astrophys. J. 471, 82 (1996)|^ivyas^-ph/9602079| 

[55] T. Sjostrand, S. Mrenna and P. Skands, JHEP 0605, 026 (2006) [; arXiv:h ep-ph/0 6031751 . 

[56] J. F. Beacom, N. F. Bell and G. Bertone, Phys. Rev. Lett. 94, 171301 (2005) 
| |arXiv:astro-ph/0409403l . 

[57] M. R. Buckley and D. Hooper. larXiv: 1004. 1644 1 Thep-phl. 

[58] R. Rando and f. t. F. Collaboration. arXiv: 0907 .06261 [astro-ph.IMl. 

[59] T. L. Collaboration. arXiv: 1003^8951 [astro-ph.CQl. 

[60] R. Aloisio, P. Blasi and A. V. Olinto, Astrophys. J. 601, 47 (2004) [ |arXiv:astro-ph/0206036| l. 

[61] A. R. Zentner, A. A. Berlind, J. S. Bullock, A. V. Kravtsov and R. H. Wechsler, Astrophys. J. 
624, 505 (2005) HarXiv:astro-ph/041 1586[ . 

[62] O. Adriani etal. [PAMELA Collaboration], Nature 458, 607 (2009) HarXiv:08 10.49951 [astro- 
ph]]. 

[63] S. W. Barwick et al. [HEAT Collaboration], Astrophys. J. 482, L191 (1997) 
[arXiv:astro-ph/9703192|; S. Coutu et ah, Prepared for 27th International Cosmic Ray Con- 
ference (ICRC 2001 ), Hamburg, Germany, 7-15 Aug 2001 

[64] M. Aguilar et al. [AMS-01 Collaboration], Phys. Lett. B 646, 145 (2007) 
narXiv:astro-ph/0703T54| . 

[65] A. A. Abdo et al. [The Fermi LAT Collaboration], Phys. Rev. Lett. 102, 181101 (2009) 
HarXiv:0905.0025l rastro-ph.HEll: D. Grasso et al. [FERMI-LAT Collaboration], Astropart. 
Phys. 32, 140 (2009) I ar Xiv:0905.0"636l rastro-ph.HEll . 

[66] J. Chang et al, Nature 456, 362 (2008). 

[67] D. Hooper, P. Blasi and P. D. Serpico, JCAP 0901, 025 (2009) HarXiv:08 10. 15271 rastro-phll : 
H. Yuksel, M. D. Kistler and T. Stanev, Phys. Rev. Lett. 103, 051 101 (2009) r arXiv:0810.2784l 
[astro-ph]]; S. Profumo. larXiv:08 12.44571 [astro-phi . 

[68] G. Kane, R. Lu and S. Watson, Phys. Lett. B 681, 151 (2009) HarXiv:0906.4765l [astro- 
ph.HE]]. 



34 



[69] D. Hooper, A. Stebbins and K. M. Zurek, Phys. Rev. D 79, 103513 (2009) [arXi v:0812.3"202l 
[hep-ph]]. 

[70] P. Brun, T. Delahaye, J. Diemand, S. Profumo and P. Salati, Phys. Rev. D 80, 035023 (2009). 

[71] P. Sandick, J. Diemand, K. Freese, D. Spolyar, and K. Zurek, in preparation. 

[72] H. Yuksel, S. Horiuchi, J. F. Beacom and S. Ando, Phys. Rev. D 76, 123506 (2007) 
HarXiv:0707.0196l [astro-ph] ] . 

[73] D. Spolyar, M. R. Buckley, K. Freese, D. Hooper and H. Murayama, arXiv: 0905 .47641 [astro- 
ph.CO]. 

[74] S. K. Mandal, M. R. Buckley, K. Freese, D. Spolyar and H. Murayama, Phys. Rev. D 81, 
043508 (2010) r arXiv:09 1 1 .5 1 881 [hep-phll . 

[75] P. Sandick, D. Spolyar, M. R. Buckley, K. Freese and D. Hooper, Phys. Rev. D 81, 083506 
(2010) HarXiv:0912.0513| [astro-ph.COll. 

[76] R. Essig, N. Sehgal, L. E. Strigari, M. Geha and J. D. Simon. larXiv: 1007.4199 [astro-ph.CO] . 



[77] G. Bertone, Phys. Rev. D 73, 103519 (2006) [arXiv:astro-ph/0603148|. 

[78] A. Achterberg et al. [IceCube Collaboration], Phys. Rev. D 75, 102001 (2007) 
[ |arXiv:astro-ph/06 1 10631 . 

[79] R. Abbasi et al. [The IceCube Collaboration], Phys. Rev. Lett. 103 (2009) 221102 
HarXiv:0911.2338i rastro-ph.HEll. 

[80] S. Toscano and f. t. A. collaboration. larXiv:0908l)864l [astro-ph.IMl . 

[81] E. Thrane et al, Astrophys. J. 704, 503 (2009) HarXiv:0907. 15941 rastro-ph.HEll . 

[82] C. Amsler et al. [Particle Data Group], Phys. Lett. B 667, 1 (2008). 

[83] C. Wiebusch and f. 1. 1. Collaboration. larXiv: 0907 .22631 [astro-ph.IMl . 

[84] R. B. Larson, |arXiv:astro-ph/9912539 

[85] V. Bromm, P. S. Coppi and R. B. Larson, Astrophys. J. 564, 23 (2002) 
JarXiv : astro-ph/0 1 02503[ ; T. Abel, G. L. Bryan and M. L. Norman, Astron. Astrophys. 
540, 39 (2000) [ arXiv:astro- ph/0002~B5]|, Science 295 (2002) 93 QarXiv: astro-ph/0 1 12088] ; 
F. Nakamura and M. Umemura, Astrophys. J. 548, 19 (2001) [arXiv:astro-ph/00 10464]; 
N. Yoshida, K. Omukai, L. Hernquist and T. Abel, Astrophys. J. 652, 6 (2006) 
narXiv:astro-ph/0606T06| . 

[86] R. Bernabei et al, Eur. Phys. J. C 67, 39 (2010) f arXiv: 1002. 10281 rastro-ph.GAll . 

[87] J. I. Collar. [arXiv:1006.203T1 rastro-ph.COl. 



35 



[88] C. Savage, G. Gelmini, P. Gondolo and K. Freese. larXiv: 1006.09721 rastro-ph.COl. 

[89] D. Hooper, J. I. Collar, J. Hall and D. McKinsev. larXiv: 1007. 10051 fhep-phl. 

[90] P. Sorensen, "Giving up S2/S1 discrimination in exchange for a lower energy threshold," 
presented at Identification of Dark Matter 2010, Montpellier, France, July 26, 2010. 



36 



